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Abstract 

We discuss in detail signals for CP violation in the Higgs boson and tau-slepton sectors 
through the production processes ^ 4 ' \i~ — > Tj~f,- , where i,j = 1,2 label the two r slepton 
mass eigenstates in the minimal supersymmetric standard model. We assume that the 
soft breaking parameters of third generation sfermions contain CP violating phases, which 
induce CP violation in the Higgs sector through quantum corrections. We classify all the 
observables for probing CP violation in the Higgs boson and r slepton sectors. These 
observables depend on the initial muon beam polarization, where we include transverse 
polarization states. If the heavy Higgs bosons can decay into tau slepton pairs, a complete 
determination of the CP properties of the neutral Higgs boson and r-slepton systems is 
possible. The interference between the Higgs boson and gauge boson contributions could 
also provide a powerful method for probing CP violation, if transversely polarized muon 
beams are available. We show in detail how to directly measure CP violation in the tau 
slepton system, under the assumption that the neutral Higgs mixing angles are determined 
through the on-shell production of the neutral Higgs bosons. 

PACS number(s): 11.30.Pb, 11.30.Er 
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1 Introduction 



CP violation was observed in the neutral kaon system [|IJ and it is strongly suggested by recent 
experiments in 5-meson decays || . In addition, CP violation constitutes one of the conditions 
for a dynamical generation of the cosmological baryon asymmetry ||. In the Standard Model 
(SM), which contains only one physical neutral Higgs boson, the only source of CP violation is 
the complex phase of the quark mixing matrix On the other hand, models with multiple 
Higgs doublets can have additional sources of CP violation, e.g. neutral and/or charged Higgs 
bosons of no definite CP quantum number. 

Supersymmetry (SUSY) is now widely regarded to be the most plausible extension of the SM; 
among other things, it stabilizes the gauge hierarchy || and allows for the Grand Unification 
of all known gauge interactions ||. Of course, supersymmetry must be (softly) broken to be 
phenomenologically viable. In general this introduces a large number of unknown parameters, 
many of which can be complex J7J. CP-violating phases associated with sfermions of the first 
and, to a lesser extent, second generation are severely constrained by bounds on the electric 
dipole moments of the electron, neutron and muon. However, it has recently been realized || 
that cancellations between different diagrams allow some combinations of these phases to be 
quite large. Even in models with universal boundary conditions for soft breaking masses at some 
very high energy scale, the relative phase between the supersymmetric higgsino mass parameter 
\x and the universal trilinear soft breaking parameter A can be 0(1) 0. If universality is 
not assumed, the phases of third generation trilinear soft breaking parameters are essentially 
unconstrained. In fact, some of these phases might be large [10] so as to provide non-SM 
sources of CP violation required for a dynamical generation of the baryon asymmetry of the 
Universe. 

Recently it has been realized that the explicit CP violation in the mass matrices of the third 
generation squarks with such possibly large CP phasesf] can induce sizable CP violation in the 



MSSM Higgs sector through loop corrections 0, [H| [14], [T5|, |16| . This induced CP violation 



in the MSSM Higgs sector can affect the phenomenology of the Higgs bosons at present and 



future colliders |L2], [13|, [Tj| [18], [19], |20], |21|j . These phases also directly affect the couplings 
of Higgs bosons to third generations sfermions. These couplings play an important role in the 
calculation of loop-induced couplings of Higgs bosons to photons and gluons as well as in the 
production of third generation sfermions at colliders |T7J . 

In the past few years a considerable amount of effort has been devoted to investigations 



of the physics potential of high energy colliders as a plausible future experimental 

program. Since muons emit far less synchrotron radiation than electrons do, a muon collider 
might be significantly smaller and cheaper than an e + e~ collider operating at the same center- 
of-mass energy y^. The main physics advantage of muon colliders compared to the conventional 
e + e~ and hadron colliders is that the larger Yukawa coupling of muons in many cases admits 



tin principle CP could also be violated in the SM by the QCD 9— term, but bounds on the electric dipole 
moment of the neutron force #qcd to be very small ^ 10~ 10 . 

Hi tan/3 is large, the third-generation CP phases are constrained by the electric dipole moments of the 
electron and neutron because in this case their contributions at the two-loop level could be sizable fj^TJ . 
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copious production of Higgs bosons as s— channel resonances, allowing to perform precision 
measurements of their properties p2|, 53, E3, IT71 |25|, |26|. In particular, they can allow us to 



search for CP-violation in the couplings of Higgs bosons to not only heavy SM fermions but 
also to third generation sfermions, which have large Yukawa couplings. 

It is difficult to probe these CP phases through processes controlled by gauge interactions, 
where large CP-odd asymmetries can emerge only if some sfermion mass eigenstates are closely 
degenerate, with mass splitting of the order of the decay width, in which case flavor or chirality 
oscillations can occur [p_7|, ^9|. On the other hand, in the MSSM CP-violating phases can 
appear at tree level in the couplings of a single sfermion species to neutral Higgs bosons. These 
phases can give rise to large CP-odd asymmetries regardless of sfermion mass splittings. Here 
we focus on f pair production. Unlike sfermions of the first two generations, f's generally 
have sizable couplings to heavy Higgs bosons even if the latter are much heavier than mz- 
Furthermore, unlike for b and t production the charge of a produced f is usually readily mea- 
surable; this is necessary for the construction of most CP-odd asymmetries. Finally, in most 
models sleptons are significantly lighter than squarks, making it easier to study them at lepton 
colliders. 

In the present work we extend the previous work |T7[ significantly in order to discuss in more 
detail signals for CP violation in the production processes — > fffj~ where i,j = 1,2 label 
the two t slepton mass eigenstates. We work in the framework of the MSSM with explicit CP 
violation. More specifically, we present a general formalism and a detailed analysis of the effects 
of the CP-violating Higgs boson mixing, and of the CP phases in the tau slepton mass matrix, 
on the polarized cross section for tau-slepton pair production. We consider both longitudinal 
and transverse polarization of the initial muon and anti-muon beams. A detailed procedure is 
suggested for measuring CP violation directly in the tau slepton system, under the assumption 
that the neutral Higgs mixing angles are determined through the on-shell production of the 
neutral Higgs bosons as demonstrated in Ref. [pl]| . The polarization observables allow for the 
complete determination of the CP properties of the neutral Higgs boson and r slepton systems. 
The interference between Higgs and gauge boson exchange contributions plays a crucial role in 
this analysis. 

The remainder of this article is organized as follows. Section |2| is devoted to a brief review 
of the mass spectra and mixing patterns of the neutral Higgs bosons and tau sleptons, focusing 



on the role of the CP phases |L3], [15], [Tq| . In Sec. y we present the helicity amplitudes of 
the production of a r-slepton pair in collisions with polarized muon beams. We also 

give a complete classification of the CP and CPT transformation properties of the polarization 
observables. In Sec. |] we show how to extract the rate and polarization asymmetries by taking 
appropriate configurations of the muon and anti-muon spins. We then perform a detailed 
numerical analysis for a representative SUSY parameter set so as to get a concrete estimate 
of the relative usefulness of those observables. Section |] is devoted to a brief summary of our 
findings and to conclusions. 
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2 Supersymmetric Particle Mixing 



Bounds on CP violating flavor changing processes tell us that the CP violating phases associated 
with flavor off-diagonal trilinear scalar interactions must be strongly suppressed ]30]. We 



therefore neglect all these flavor changing CP violating phases in the present work, so that 
the scalar soft mass matrices and trilinear parameters are flavor diagonal and the complex 
trilinear terms are proportional to the corresponding fermion Yukawa couplings. Clearly the 
Yukawa interactions of the top and bottom (s) quarks play the most significant role in radiative 
corrections to the Higgs sector. Furthermore, for large values of tan j3 the fi and tr states are 
expected to mix strongly. 



2.1 CP— violating neutral Higgs boson mixing 



In this section, we give a brief review of the calculation |T3) of the Higgs boson mass matrix 
based on the full one-loop effective potential, valid for all values of the relevant third-generation 
squark soft-breaking parameters. The CP-violating phases in the top and bottom squark 
sectors cause scalar-pseudoscalar mixing at one-loop level. 

The MSSM contains two Higgs doublets Hi, H 2 , with hypercharges Y(H\) = —Y(H 2 ) = 
— 1/2. Here we are only interested in the neutral components, which we write as 

#?=-^(0i+ifli); # 2 ° = ^(0 2 + m 2 ), (1) 

where <pi t2 and aia are real fields. The constant phase £ can be set to zero at tree level, but 
will in general become non-zero once loop corrections are included. 

The mass matrix of the neutral Higgs bosons can be computed from the effective potential 







Vffiggs = \m\ {4>\ + a\) + \m\ + a 2 2 ) - \m{ 2 \ (0i0 2 - aia 2 ) cos(£ + 6 U ) 

+ \m\ 2 \ {<Pia 2 + <p 2 ai) sin(£ + 9 U ) + f V 2 + ^Str [M A (log ^ - §)] . (2) 

We have allowed the soft breaking parameter m\ 2 = \m\ 2 \e iei2 to be complex, and we have 
introduced the quantities 

T-> / 2 i 2 ±2 2 -2 9 2 + 9 2 f \ 

V = (j) 2 + a 2 -a x ; g = , (3) 

where the symbols g and g' stand for the SXJ(2) L and U(l)y gauge couplings, respectively. Q 
in Eq. @ is the renormalization scale; the parameters of the tree-level potential, in particular 
the mass parameters m\, m 2 and m\ 2 , are running parameters, taken at scale Q. The potential 
(H) is then independent of Q, up to two-loop corrections. 

The matrix M. in Eq. (0) is the field-dependent mass matrix of all modes that couple to 
the Higgs bosons. The by far dominant contributions come from the third generation quarks 
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and squarks. The (real) masses of the former are given by 



m. 



^\h h \ 2 



(<t>l + al) ; 



m. 



(02 + 4) 



(4) 



where h and h t are the bottom and top Yukawa couplings. The corresponding squark mass 
matrices can be written as 



Ml 



( 2 

1 m~ 



+ m 2 -\{g 2 -^)V -h* t [A* t {Hiy+nHl 



-h t 



A t H$ + ii* (iZ?) 



m~ + mf 



Ml 



V 



-h* 



\ 



-h 



AbHl + /i* (H° 2 ) 



m 2 ~ + m 2 + alv 



(5) 



Here, H® and H® are given by Eq. ([I]) while m\ and m 2 are as in Eq. (PS) and T> has been 
defined in Eq. (Q). In Eq. ([5]) m~, m~ and m~ are real soft breaking parameters, A and A t 

are complex soft breaking parameters, and \x is the complex supersymmetric Higgs(ino) mass 
parameter. 

The mass matrix of the neutral Higgs bosons can now be computed from the matrix of 
second derivatives of the potential @, where mf, m\ and m\ 2 sin(£ + # 12 ) are determined 
by the stationarity conditions^ dV/dipi = dV/d(p2 = dV/dai = 0. The massless state G° = 
ai cos (3 — a 2 sin (3 is the would-be Goldstone mode "eaten" by the longitudinal Z boson. We 
are thus left with a squared mass matrix M 2 H for the three states a = a± sin/5 + c^cos/?, 0i 
and 02- This matrix is real and symmetric, i.e. it has 6 independent entries. CP violation is 
caused by the appearance of nonvanishing a<$>\ or a<p2 entries of this mass matrix. The sizes of 
these entries are controlled by the quantities 



A; 



Sm(y4 i /ie^) 

2 2 

mf — mf 

*2 tl 



Qm(A b fie^) 

2 2 

mf — mf 

02 01 



(6) 



Explicit expressions for M\ can be found in Ref. [15||. 

Since M\ is real and symmetric, it can be diagonalized with a three-dimensional orthogonal 
rotation matrix 0\ 



1 a \ 




( Hi \ 


01 


= o 




V 02 J 




KH 3 ) 



(7) 



Our convention for the three mass eigenvalues is mut < ^h 2 — mfl 3 . The loop-corrected 
neutral Higgs boson sector is thus determined by fixing the values of the following parameters: 
rriA, which becomes the mass of the CP-odd Higgs boson if CP is conserved, and tan/3 fix the 



'''The condition dV/da 2 = gives the same constraint as dV/dai = JO 
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tree-level Higgs potential; and //, A t} Af, and the soft-breaking third generation sfermion masses 
tuq, my, and m^, which fix the third generation squark mass matrices. After minimization of 
the potential the rephasing invariant sum 12 + £ of the radiatively induced phase £ and the 
phase 9 12 of the soft breaking parameter m\ 2 is no longer an independent parameter. [] The 
physically meaningful CP phases in the Higgs sector are thus the phases of the re-phasing 
invariant combinations Atfie^ and A^fie^ appearing in Eqs. (^). We refer to Ref. [13] for 



further details on neutral Higgs boson mixing in the presence of CP-violation. This mixing 
changes all the couplings of the Higgs fields so that the effects of CP violation in the Higgs 



sector can be probed through various processes |Tq, 19, 20, 21 



2.2 Tau slept on mixing 

The mass matrix squared M. T for the left-/right-handed tau sleptons is given by 




M * = U : - 1 - m 

where the matrix elements are defined as 

X T = m\ +ml + \{m\ - 2m^) cos 2(3 , 

O O/O O \ „ 

Y T = m f + m T + (m w — m z ) cos2p , 
Z T = m T | A* + ji e^ tan (3 \ , 

4> T = arg(-A T - ix* e~ i€ tan/3) , (9) 

where m? are the left /right-handed soft-SUSY-breaking tau slepton masses squared, respec- 
tively. 

The mass matrix squared Ai? can be diagonalized by a unitary transformation U T : 

UlMlU T = dmg(m 2 fi ,ml), (10) 

with the mass ordering m fl < convention. The diagonalization matrix can be param- 

eterized as 





Ur = . n „ , (11) 



*As discussed in flTSl , ^ and 9\2 are not separately physical parameters. For example, one or the other can 
be set to zero in certain phase conventions for the fields. Similar remarks hold for the phases of A t , Af, and \x. 
Altogther there are only three rephasing invariant (i.e. physical) phases, which we write as 612 +£, arg(A t /^e^) 
and arg(^4f,/ie J £). The minimization of the potential fixes one of these combinations, leaving two independent 
physical phases as free input parameters. 
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taking the ranges, — n/2 < 9 T < and < T < 27r. The tau slepton mass eigenvalues and 
mixing angles are then given by 

sin2fl T = -2 4 1 , cos2fl T = - Xt ~ Yt . (12) 

2 

The phenomenological parameters M f and A T are related to the tau slepton masses as 

¥ 2 = rr& + m l _ X T + Y T 

f ~ 2 2 ; 

Ar = m 2 f2 - ml = ^(X T - Y T y + 4Z T 2 . (13) 
Clearly tau slepton left-right mixing is enhanced for large tan/3 and large 



3 Tau Slepton Pair Production 

In this section, we present all the Feynman rules needed for the production process — > 
f~fj~ = 1,2), derive the production amplitudes, and describe how to obtain the polar- 
ized cross sections with initial muon and anti-muon polarizations. We then classify all the 
polarization and rate asymmetries according to their CP and CPT transformation properties, 
and discuss the CPT-even and -odd combinations of the neutral Higgs boson propagators, 
neglecting the Z boson width. 



3.1 Feynman rules 

The couplings of the 7 and the neutral gauge boson Z to fermions in the MSSM is described 
by the same Lagrangian as in the SM: 

C vff = -eQ f h,f A" - -2-/ 7|t [ (T /3 - Q f s 2 w )P^ - Q f s 2 w P + ] f Z» , (14) 

Cw L J 

with e = gsw and the chirality projection operators P± = (1 ± 7s)/2. Qf is the electric charge 
of fermion / in units of the proton charge. The couplings of the neutral Higgs bosons to leptons 
and quarks are described by the Lagrangian 

h 3 hd 3 

£aff = 7= ? [ °2fc - is /3 0^75] £H k -= Y d [ 2k - isp lk ^ 5 ] dH k 

v 2 k=1 V2 k=1 



K 3 

7= u [ 3k - icp Oifc7 5 ] uH k . 

V2 k=1 



(15) 
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Here hi, hd, and h u are the lepton and quark Yukawa couplings: 

gnu gm d gm u 
h i = —fK , h d = -j= , h u = — = , (16) 

respectively, with cp = cos (5 and sp = sin (5. It is then clear that all the neutral Higgs bosons 
couple dominantly to the third generation fermions t, b and r and they couple to a muon about 
200 times more strongly than to an electron - the primary reason for having a muon collider. 

The couplings of the 7 and the neutral gauge boson Z to tau sleptons in the MSSM are 
described in the mass eigenstate basis by the Lagrangian 

Cvnf, = * eA^j (f * X fj) -i — Z" Q% {f* X fj) , (17) 

Cw 

where Qfj are expressed in terms of the tau slepton mixing matrix U T as 

<% = *w8ii-\(Ur)\&J T )v. (18) 

On the other hand, the Feynman rules for the Higgs boson couplings to tau slepton pairs, 
involving all the mixing and phases and including the phase rotations of the scalar tau fields, 
can be written as 

Chtt = — - — - O a k{U T )*p i (U T ) jj V a ;f3.y H k f*fj = — - — - Vk-ij H k ?*fj (19) 

Cw Cw 

where a labels the three neutral Higgs boson interaction eigenstates {a, 0i,02j-, and {/3,7} 
denote the chiralities {L, R} of the f interaction eigenstates. The chiral couplings V a;/ 3 7 for the 
scalar tau leptons can be obtained in a rather tedious but straightforward way as 



V, "' 



i(A* tan f3 - ne^) 



2m| y -HArttmP- n*e-*) 



2 _ 1 1 m T mT . A Z. 
W 2 ' m 2 z cj 2m z c j 



m T A T 2 I rK 

2^PI S W+m£cT 



Vfo = -Sp 



/ ,.2 _ 1 m T /ie'5 
" W 2 Zrn^cpsp 

in T ii*e-'''- _„2 



(20) 



These 2x2 matrices determine the dimensionless couplings Vk-ij through the relation: 

Vk-ij = O ak (UrYpt (U T ) yj Va-fa. (21) 
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Figure 1: The mechanisms contributing to the process — > fffj : three spin-0 neutral- 
Higgs boson exchanges and spin-1 7 and Z exchanges. The indices i,j are 1 or 2, and the 
index k is 1,2 or 3. 



3.2 Production amplitudes 

As shown in Fig. 1, the matrix element Aiij for — > fffj~ receives contributions from 7 
and Z exchange as well as from the exchange of the neutral Higgs bosons of the MSSM. It can 
be written as: 



Mij = — J2 \ Z *ij [ v (p)(^i- tj) P Mp)) + m w H aij [v(p) P a u(p)} X , (22) 

S ,v=4- I. > 



-,2 

a=± 

where a — ±, P± — (1 ± 7s)/2, and the vector and scalar chiral couplings are given by 



Z-ij = Sij + D z (s) Sw y 2 Qfj , 
c w b w 

Z +ij = 5 i:j + D z (s)— Q§ , 
c w 

H-ij = —Y,i 2 k ^2 ' [^ 2fc + ispOik] , 
c w s w 

H +ij = -Y„ ^fi- V k;tJ [0 2k - ispO lk ] . (23) 
c w s w 

The coupling Y" M , and the Z and H k propagators D z and D Hk (k = 1, 2, 3) are: 

Y = m » 

2m w cp ' 

D z(s) = * r , 

s — m z + IVTLzV z 
Dh ^ S) = s m* 1 4m V " (24) 



9 



Denning the polar angle of the flight direction of the tau slepton ff with respect to the fi~ 
beam direction by (see Fig. 2), the explicit form of the production amplitude (22) can be 
evaluated in the helicity basis by the 2-component spinor technique of Ref. J32JP]. We denote 
the \T~ and /i + helicities by a and a, with a = + and — standing for right- and left-handed 
particles, respectively. Neglecting the muon mass in the kinematics, the helicity amplitudes 



Mij(ar,a) 



(aa) 



i.i ■ 



(25) 



read 



m w 



H 



's 
m w 



H 



("+)« 



•Z+ij (3 sin 9 , 
-Z_ij (3 sin 6 . 



(26) 



Here (3 = A 1 / 2 (l, m?_/s, m?+/s), with X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + yz + zx), describes 

the 3-momentum of the scalar r leptons in the center-of-mass frame. Note that in the limit of 
negligible kinematic mass of the muons vector boson exchange only contributes to configurations 
where \i~ and /i + have opposite helicities, whereas Higgs boson exchange only contributes to 
configurations with equal fi + and fi~ helicities. This implies that Higgs and gauge boson 
exchange contributions can only interfere with each other if at least one of the muons has 
nonvanishing transverse polarization, which can be understood as coherent superposition of 
left- and right-handed helicity states |32| ; see Eqs. fl3lf ) and (|||) below. 




Figure 2: Schematic depiction of the production plane with the polar scattering angle G. The 
transverse polarization vectors Pt and Pt have azimuthal angles a and a with respect to the 
event plane, respectively. 



Before turning to the calculation of the cross section, let us briefly describe the symmetry 
properties of the helicity amplitudes. The CP transformation leads to a relation among the 
transition helicity amplitudes: 

(vvjij " (-*,-*)„; (27) 

^Our convention for the momentum-dependent Weyl spinor for fermions going in —z direction differs by an 
overall sign from that of Ref. Q 
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note that a — a gives the total spin, if the fi momentum direction is used as quantization axis. 
Equivalently one has for the scalar and vector helicity amplitudes: 

(±±)« " -<=FF>ii, (±T)ij " +(±T)ii- (28) 

Eqs. (^8|) , ( ^3|) and (|26|) show that CP violation in diagonal channels (i = j) can only occur in 
the presence of nonvanishing pseudoscalar couplings, since the contributions oc 0^ in Eq. ( |23| ) 
are the only ones that violate relation ( p8|) . We will see in the next section that observable CP- 
odd effects can only arise in diagonal f pair production if a pseudoscalar coupling interferes with 
one of the other couplings. On the other hand, all four sets of couplings (scalar, pseudoscalar, 
vector, and axial vector) can give rise to CP violation in non-diagonal f pair production (i =f j). 

Another useful classification is provided by the so-called "naive" time reversal T. Like 
the proper T-transformation, this transformation changes the directions of all 3-momenta and 
spin vectors, i.e. it leaves helicities unchanged. However, unlike a T transformation, a T 
transformation does not interchange initial and final states. In our case a T transformation 
simply corresponds to taking the complex conjugate of the helicity amplitudes. Under the 
combined CPT transformation the helicity amplitudes therefore transform as follows: 

(±±>« " -<=F=F>i, (±=F>« " +<±=F>£. (29) 

We note that tree-level CPT violation is only possible in the presence of finite Z or Higgs boson 
widths. It is therefore very useful to analyze the CP and CPT properties of any physical ob- 
servable simultaneously, so as to investigate not only CP violation itself but also its dependence 
on the real or imaginary part of some product (s) of propagators. 



3.3 Polarized production cross sections 

The production cross section for general (longitudinal or transverse) beam polarization can be 
computed either using standard trace techniques (employing general spin projection operators), 
or from the helicity amplitudes by a suitable rotation from the helicity basis to a general 
spin basis. In the former case, neglecting the muon mass in the spin projection operators we 
can obtain the following approximated form for the fi T projection operators: 

-(p + m)(l +75 f) — > ~(1 + PLlb) P + -75-Pr(cosa fa + sina fa) , 
-($-m)(l + 7 5 f) — ► -(1 -Pz/y 5 ) P+ ^ 5 P T (cosa fa + sina fa) p. (30) 

Here s M is the spin 4-vector, and rii^ and rii^ are purely spatial vectors satisfying rii ■ rij = —5ij 
and rii ■ p = rii ■ p = 0. Pl and Pl are the longitudinal polarizations of the \T and fi + beams, 
while Pt and Pt are the degrees of transverse polarization with a and a being the azimuthal 
angles between the transverse polarization vectors and the momentum vector of f~ as shown 
in Fig. 2. We also note that P| + P| < 1 and P| + P| < 1. 
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Equivalently the polarization weighted matrix element squared in the helicity basis is given 



by 



Y,\M\ 2 = £ M a - a Ml, &l p- a ,pt- a , = Tt M P + M^p- T 



(31) 



where M. a5 (a, a — ±) denotes the helicity amplitude for any given production process p (a)p + (a) 
X and the 2x2 matrices p T are the polarization density matrices for the initial p T beams: 



1 1 + P L P T e- i0 

2 \ P T e ia 1 -Pt 



1+P 



-P T e 



l ~P T e ta 

ia i _ p 



(32) 



Applying the projection operators (|30|) or/and evaluating the trace (|31| ) leads to the follow- 
ing matrix element squared for polarized p + p~ — > fffj" production: 



1 - P L P L ) C x [ij] + [P L - P L ) C 2 [ij] 
+ (l + P L P L ) C 3 [tj] + (P L + P L ) C 4 [tj] 
+ (Pp cos a + Pt cos a) C 5 [ij] + (Pr sin a + Pr sin a) C 6 [ij] 
+ (Pr cos a — Pt cos a) CV[ij] + (Pr sin a — Pt sin a) Cg[ij] 
+ (P L Pr cos a + P l Pt cos a) C 9 [ij] + (P L Pr sin a + P l Pt sin a) C w [ 
+ (PlPt cos a — PlPt cos a)Cn [ij] + (PlPt sin a — PlPt sin a) Ci 2 [ 
+ PrPr [cos(a + a) C\z[ij] + sin(a + a) Cu[ij\\ 
+ PtPt [cos(a — a) C^[ij] + sin(a — a) C\q[ 



The coefficients C n (n — 1 - 16) are defined in terms of the helicity amplitudes by 



(33) 



Ci 

c 3 
c 5 

cv 
c 8 

C 9 



K+-)I 2 + I(-+)P 
K++)l 2 + K— )P 



G 2 



l(+-)l 2 -K-+)l 2 
l(++)l 2 -K— >l 2 



He ((++)_(__))((_+)_ (+_))* 



9ta ((++>- (—»((-+> + (+-»• 

«e «++> + <— »«-+> + <+-»*, 
Sfen «++> + (— »((-+>-(+-»• 

»e«++> + (—))((-+) -<+-»*, 
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Cio 
Cn 

Cl2 
C13 
C15 



J9ta ((—} + (++))((-+} + (+-})*, 
ifce «—)-<++»«-+) + <+-))*, 

i9(m ((++>- (—»((+->-(-+»*, 



--«e [(-+)(+- 



C 14 = -3m [<-+)(+-)*], 



i«e [(—>(++>*], Ci6 = ~ 3m [<—><++>*]. 



(34) 



The production cross section is then given in terms of the distribution E,,- in Eq. 



by 



d. 



a 



d cos O rf$ 



T- T ■ 



4 s ^ « ' 



(35) 



with /3 as in Eqs. (p6|). The dependence of the distribution Ejj on the azimuthal angle $ of 
the production plane is encoded in the angles a and a. If the azimuthal angle $ is measured 
with respect to the direction of the pT transverse polarization vector, the $ dependence can 
be exhibited explicitly by taking for the angles a and a 



a 



a 



7] - $ 



(36) 



where rj is the rotational invariant difference a — a of the azimuthal angles of the fi + and \l~ 
transverse polarization vectors. The polarization coefficients can be classified according to their 
correlation patterns of the vector and scalar contributions as follows: 



Vector and vector (VV) correlations 
Scalar and scalar ( SS ) correlations 
Scalar and vector ( SV) correlations 



C3, C*4, Cl5; Cl6; 

5; ^6; Wj ^9) WO) WI5 W2- 



(37) 



Among the sixteen observables, only six observables - the 2 VV observables {C\, C2} and the 
4 SS observables {C3, C4, C15, C\q\ - can be measured independently of the azimuthal angle 
but the other 10 observables require the reconstruction of the production plane. This is true in 
particular for the observables involving scalar-vector correlations; as already stated in Sec. 3.2, 
these observables are only nonzero if at least one of the incoming beams has nonvanishing 
transverse polarization. 



3.4 Observables with definite CP and CPT properties 

The CP and CPT transformations act on the polarization vectors of the initial muon beams 
according to the simultaneous exchanges: 



Pr 



CP 



CPT 



-Pl, 
-Pl, 



Pr 



CP 



Pl 



P T « P T : 



a 



a 



CP 



a: 



CPT 



-> —a . 



(38) 
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The CP and CPT relations (^) of the polarization vectors lead to the following classification 
of the terms in Eq. (p3|): 7 CP-even and CPT-even, 3 CP-even and CPT-odd, 3 CP-odd and 
CPT-even, and 3 CP-odd and CPT-odd polarization factors, as shown in Table |]. 

Table 1: The classification of the correlations of the muon and anti-muon polarization vectors 
that appear in Eq. ftS^j ) according to their CP and CPT properties, based on the CP and CPT 
relations (p^j. 



CP 


CPT 


Polarization factors 


even 


even 


1-P L P L , Pl-Pl, 1 + PlPl, Pt cos a + P T cos a 
PlPt cos a — PlPt cos a, P T P T cos(a + a), P T P T cos(a — a) 


odd 


Pt sin a + Pt sin a, PlPt sin a — PlPt sin a, P t ?t sin(a + a) 


odd 


even 


Pt sin a — Pt sin a, PlPt sin a + PlPt sin a, PtPt sin(a — a) 


odd 


Pl + Pl, Pt cos a — Pt cos a, Pl-Pt cos a + Pl-Pt cos a 



The coefficients C n [y] (n = 1-16) corresponding to these polarization factors are bilinears 
in Z±ij and H±ij. For later convenience we define the asymmetric combinations [C n [12]] and 
the symmetric combinations {C n [ij']} (ij = 11, 12,22) as 

[C n [12]] = 1 (C„[12] - C n [21] ) , {C n [ij]} = i (C„[y] + C n [jz] ) . (39) 

Clearly the coefficients multiplying CP-even polarization factors (the first group in Table 1) can 
give rise to CP violation only through the antisymmetric combinations [C n [12]], which describe 
rate asymmetries. In contrast, all symmetric combinations {C n [r/]} of the coefficients of the 
second group in Table 1, which have CP-odd polarization factors, can contribute to CP-odd 
polarization or azimuthal angle asymmetries; these can be probed for three different CP-even 
final states (f^f^, i — 1, 2 and the sum of rffjj 1 " and f+f^ production). Based on these 
observations one can classify all the 64 distributions according to their CP and CPT parities 
as shown in Table We will show in Sec. 4.2 that all 64 coefficients appearing in Eq. (|33D can 
be extracted independently, if the polarization of both beams can be be controlled completely 
and the final tau slepton pair is identified. 

We thus see that one can define 10 CP-violating rate asymmetries involving anti-symmetric 
combinations in the f indices and 6 polarization/angle asymmetries involving symmetric com- 
binations in the f indices; we recall that the latter can be studied for three different final states, 
leading to a total of 28 different asymmetries. More detailed explanations of the definition of 
the rate and polarization asymmetries and of the extraction of the asymmetries by adjusting 
the muon and anti-muon polarizations will be given in the next section. 
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Table 2: The classification of the coefficients C n [ij] (n = 1-16, ij = 11, 12, 21, 22 ) according 
to their CP and CPT properties. The underlined observables can be constructed without direct 
reconstruction of the scattering plane. Here, {C n [ij}} = (C n [ij] + C n [ji])/2 and [C n [12]] = 
(C n [12]-C n [21])/2. 



CP 


CPT 


Coefficients (ij = 11,12,22) 


even 


even 


{c x [ij\}, {c 2 [ij}}, {c 3 [ij}}, {c 5 [ij}}, {CiiM}, {c 13 m 

{C l5 [ij}}, [C 4 [12]}, [C 7 [12]], [C 9 [12]] 




odd 


{C 6 [ij\}, {C 12 [ij\}, {C u [ij\}, [C 8 [12]], [C 10 [12]], [C 16 [12]] 


odd 


even 


{C 8 [ij}}, {C 10 [ij}}, {C 16 [ij}}, [C 6 [12]], [C 12 [12]], [C 14 [12]] 




odd 


{C^ij}}, {Cr[ij}}, {C 9 \ij]h \Cmi [C a [12]], [C 3 [12]] 
[C 5 [12]], [C u [12]], [C 13 [12]}, [C 15 [12]] 



Any CPT-odd observable requires some CP-preserving phase. At tree-level this can only 
be provided by the Z and Higgs boson widths in the tau slepton pair production processes. 
Therefore, if the Z width is neglected, in which case Z±^ is Hermitian, the following six 
quantities involving only the 7 or Z boson couplings to muons and tau sleptons vanish : 

[d[12]] = [C 2 [12]] = [C 13 [12]] = C 14 [11] =C 14 [22] = 0. (40) 

Furthermore, Eqs. (|18D , 

[C7 14 [12]] = {C 14 [12]} = 0, (41) 

even for nonvanishing V z- All other CPT-odd quantities can get contributions from some Higgs 
boson width, as will be clarified in the next section. 

3.5 CPT— even and odd combinations of Higgs boson propagators 

The behavior of each observable under a CPT transformation plays a crucial role in determining 
the interference pattern of the Higgs boson contributions, among themselves and with the 7 
and Z boson contributions to this observable. Every CPT-even (odd) observable involving only 
scalar contributions depends on the real (imaginary) parts of the products D Hk D* Hi of Higgs 
boson propagators; note that only k 7^ I gives a nonzero imaginary part. We already showed 
that the observable Ci 4 involving only the vector contributions is negligible, see Eqs. ( fjQ|) and 
(PI) . In the approximation V z = all the vector boson contributions Z_ and Z + in Eq. ( Effi) are 



23) and (pi) show that 
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hermitian. As a result, the interference terms between scalar and vector contributions depend 
on the real or imaginary parts of the Higgs boson propagators Dn k themselves. 

This discussion shows that it is convenient to introduce abbreviated notations for the real 
and imaginary parts: 



Shi = Me[D Hk D*\ = s 



2 (s-rn% k )(s-m 2 Hi )+m Hk m Hl T Hk T Hl 



o ( s — m 1r ) m uX 'h, — {s — mir) mjfj ' H , , x 

It is worthwhile to note two points: (a) the denominator reveals a typical two-pole structure so 
that the Higgs boson contributions are greatly enhanced near the poles; (b) the numerator of 
Ski is negative between two resonances with a mass splitting larger than their typical widths but 
positive otherwise. In contrast, for realistic widths of the MSSM Higgs bosons the numerator of 
T) k i is always positive (negative) if >mu k (jnH k > m Hi)- Moreover, it increases (decreases) 
linearly with s if tuh^Hi > (<) m H k ^H k - 

As noted earlier, in the approximation V z = the SV interference terms depend on the real 
and imaginary parts of each Higgs boson propagators D Hk , which are given by 

**[Dn k ] = 7 ^fegN r2 , ^m^J = S ^fj\ r2 , (43) 

(s-m 2 Hk y + m 2 Hk r 2 Hk (s-m 2 Hk ) 2 + m 2 Hk r 2 Hk 

respectively. Note that the real parts change their sign whenever the cm. energy crosses the 
corresponding pole, but the imaginary parts are always negative. 

Combining the coefficients from the mixing matrix elements with the corresponding propagator- 
dependent factors enables us to straightforwardly understand the qualitative y^i dependence 
of the combinations of coefficients C n [r/'] listed in Table 2. However, the energy dependence of 
the corresponding properly normalized asymmetries is also affected by the energy dependence 
of the normalization factors, which can be nontrivial. We will discuss these issues further in 
Sec. 5, where we present numerical results. 



4 Polarization and rate asymmetries 

In this section, we present the explicit forms of all the VV, SS, and SV polarization and rate 
correlations and then explain how to extract the corresponding asymmetries by adjusting initial 
beam polarization and, in some cases, inserting weight functions of the polar angle in the phase 
space integral. For completeness we include the CP-even correlations listed in the upper half 
of Table 2. 



4.1 Expressions for the correlations 
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4.1.1 Vector-vector correlations 

The VV correlations are very useful for reconstructing the tau slepton system independent of 
the Higgs boson system. In principle there exist 9 non-trivial CP-even and CPT-even VV 
observables ({Ci[ij]}, {C 2 [«j]}, {Ci 3 [ij]}) and 1 CP-odd and CPT-even observable [Ci 4 [12]]. 
However, we already saw that this latter observable vanishes. Since the CP-violating phase 
<fr T in the tau slepton system affects neither the masses of the f eigenstates nor the absolute 
values of their gauge couplings, it cannot be directly measured through VV correlations only.^ 
However, the absolute values of X T , Y T and Z T defined in Eq. (^) can be determined from the 
remaining 9 CP-even and CPT-even observables, together with the measurement of the two 
tau slepton masses. As a result, if tan/3 is determined (or constrained to be large) from other 
processes (e.g., chargino pair production in e + e~ collisions the above determination can be 
used to obtain the left/right soft-breaking tau slepton masses m fL R . The absolute magnitude 
of \A* + /ie*£ tan/3| can be determined even if tan/3 is unknown. However, it is necessary^ to 
consider different observables to determine the phase angle <p T = arg(— A T — (i* e~^). 



4.1.2 Scalar-scalar correlations 

The scalar-scalar correlations involve not only tau slepton mixing but also neutral Higgs boson 
mixing. However, it is known that the Higgs boson mixing angles can be completely determined 
through on-shell production of each Higgs boson with polarized muons In this section we 



therefore check the possibility of reconstructing the tau slepton system through these scalar- 
scalar correlations. There are 16 such correlations: 7 CP-even and CPT-even, 1 CP-even 
and CPT-odd, 3 CP-odd and CPT-even and 5 CP-odd and CPT-odd observables. It is 
straightforward to obtain the explicit forms of the these 16 correlations: 




s kl 3fe (u :;/ r,;,..) (o 2k o 2l + s} o xk o x 



w V kl 3m (V fc; i2v;* 12 ) (o 2k o 2l + sj o lk o u 



V kl Ke (V k;ij V^) s p (0 lk 2l - 2k O u ) , 
S k i Qm iy k ;i2^*i 2 ) sp {0 lk 2i - 2k On) . 



Ski [V k;ij V^A (0 2k 2l - si lk O 



^The phase &a t does affect the f masses. However, it could be extracted from measurements of f masses 
and 6 T only if \A T \ and tan/3 are already known. In case of \A T \ at least this seems quite unlikely. 

'•The absolute value of Z T does depend on the relative phase between A* and fie 1 *. However, observables 
that only depend on \Z T \ clearly cannot determine this phase, unless \A T \ and |/x tan /3| are already known. 
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[Cis[12]] = - (^%rj ^V kl Zm(y k . 12 VZ 12 ) (0 2k 2l - s}Oi k O u ) , 

{c le [ij}} = - (jAr) ^ 5 « (\, : , ( \;:, ; ) Sf3 (o lk o 2l + o 2k o u ) , 

\ c w s wJ ls 

( Y \ 2 m 2 / \ 

[C 16 [12]] = + ~^ V ki Sm (v^Vfa) s p (0 lk 2l + 2k O u ) , (44) 

\ c w s wJ Zs 

where summation over the repeating indices k and I is understood. In deriving Eqs. fl44|) 
we have exploited the relation Vfc ; jj = (V k .ji)*. Note that any rate asymmetry, defined to be 
proportional to the difference between the [12] and [21] modes, is determined by the combination 
3m(Vfc ; i2VJ* 12 ), which is antisymmetric with respect to the Higgs boson indices k, I. In addition, 
in the MSSM the lightest Higgs boson has a mass of less than 130 GeV and a tiny width, so that 
its contribution to the tau slepton pair production is negligible. As a result, only the propagator 
combinations S 22 , S 23 , S 33 and V 23 are important in determining the y/s dependence of these 
rate asymmetries. 



4.1.3 Scalar— vector correlations 

All the SV correlations are proportional to the real or imaginary part of a Higgs propagator, so 
that they also are sizable when the cm. energy is very close to one of the heavy Higgs boson 
resonances. § 

{C 5 [ij}} = + .7^4 OnDz MDhA He (Vw Qf*) sin e, 
4V sc w s w 

[C 5 [12]] = - A m ^ A 2k D z 9m[Z) g J Sm (v k , l2 Qf*) sin6, 
4 v sc w s w 

{C e M} = - 2k $m[D Hk ] Ke (v k , tJ 5 tJ + Dz ^'^ V^QfA sin 6, 

v sc w s w \ c w s w J 

[C 6 [12]] = - ™?f "f (s 2 w - l/4)0 2k D z Ue[D Hk ] Sm (V k . l2 Qf*) sin 6, 

V sc w s W 



{C 7 [ij}} = ~ m r 2 Y ^ 2 H lk %m[D Hk } Ke ( \ ),..,,<) ; , + d / w \^ V k ^ Qf* ) sin 6, 
V sc w s w \ c w s w 



[C 7 [12]] = - (4 - l/4)sf,0 lk D z ?fte[D Hk ] Sm (V fc;12 Qf 2 *) sin 9, 

{Cs{ij}} = - A m r^ S P Ov.Dz ®e[D Hk ] He (v k;ij Qf*) sin 6, 
sc w s w 

[ Cs[12]] = + A m ^ /3 4 sp lk D z $m[D Hk ] Sm (v k;12 Qf*) sin 6, 



§$le(DH k ) = for y/s = mn k ■ However, since the mass splitting between the two heavy Higgs states is small, 
the contribution from the second heavy Higgs boson is usually still sizable at this point. 
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{CM} = + A m r^ s ? °^ Dz <M D B h ] 3«e (V k;ij Qf*) sin 6, 
4 V sc w s w 

[ C 9 [12}} = + A m ™^ s p lk D z Me[D Hk ) 3m (V k;12 Qf*) sin 9, 

{C w m = + ^4 S P lk ne[D Hh ] Ke [V^j + d / w ~ X J\ ^ Qf*) sin 9, 
V SC W S W \ c w s w ) 

[C 10 [12\] = - (s 2 w - l/4)sp0 lk D z Zm[D Hk ] 3m (V k . 12 Qf*) sin 9, 

v sc w s w 

{CnM) = + ^4 02kMe[D Hk ] 3tfe (v^Sy + Dz^f^V^j Q?A sin 9, 
v sc ty s w V c w s w / 

[C u [12]} = - (s 2 w - l/A)0 2k D z %m[D Hk ] 3m (v k , 12 QfA sin 9, 

v sc w s w 

{C 12 [ij}} = - A m r Y / 4 2k D z Zm[D Hk ] Ke (V m Qf*) sin 9, 

[C u [12]] = - . m ™ Y /{ 2k D z Ue[D Hk ] 3m (V M2 Q§*) sin9. (45) 



The CP and CPT properties of the observables in Eq. ( |45|) can be identified by noting the 
following aspects: (i) observables involving $t,e[D Hk ] (SmjD^J) are CPT-even (CPT-odd). (ii) 
The matrix elements 2k and 0\ k are the CP-even and CP-odd component of the neutral Higgs 
boson H k , respectively, (iii) 5te (V k;i j 5ij) and 9fte (v k;i j Qff) are CP-even, while 3m (v k -i 2 Qf 2 ) 

is CP-odd. For instance, {C 5 [r/]} is CP-even and CPT-even, while [C 6 [12]] is CP-odd and 
CPT-even. Of course, this determination of the CP and CPT properties agrees with the results 
of Table 2. 



4.2 Extracting the polarization and rate asymmetries 

In this section we investigate efficient procedures for extracting the relevant SUSY parameters 
from the polarized cross sections. In actual experiments, a careful analysis of statistical and 
systematic uncertainties will be required when extracting the rate and polarization asymmetries. 
We do not attempt to estimate realistic uncertainties based on event simulation or parameter 
fitting. Instead, we are interested in the question which of the numerous quantities listed in 
the previous Subsection are most sensitive to the fundamental parameters of the theory, our 
main focus being on CP-odd phases in the interaction Lagrangian. 

The conceptually simplest asymmetry is the total unpolarized rate asymmetry, defined by 

.tot _ Jrfcos9([C 1 [12]] + [C 3 [12]]) 

R Jrfcos9({C 1 [12]} + {C 3 [12]})- 1 ' 

Unfortunately this asymmetry is usually very small even if CP is violated in both the Higgs 
and f sector. We saw in Eq. (^) that [Cx[12]] vanishes in the limit T z — > 0. Moreover, 
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Eq. (fH]) shows that [C/3[12]] is proportional to 2 k0 2 i + sin 2 fiOikOu. The contribution from 
Higgs exchange can only compete with the gauge boson exchange contributions if ^/s ~ mHi ggs - 
Given LEP search limits for fx, real f pair production will only be possible at energies s ^> m%.* 
These two requirements together imply that we need m 2 A 3> m 2 z . In this "decoupling limit" the 
lightest neutral Higgs boson eigenstate is given by Hi ~ cos/30i + sin/?0 2 ; the couplings of this 
state are very similar to that of the Higgs boson of the Standard Model. The rotation matrix 
O of Eq. fl?D is then approximately given by: 



O 



I cosa# sina# 

cos (3 sin ot H sin (3 — cos och sin (3 | , (47) 
^ sin (3 — sin a# cos (3 cos an cos (3 



where sm(2a#) 7^ signals CP-violation in the Higgs sector. The value of ajj depends in 
a complicated manner on the parameters appearing in the squark mass matrix, and on tua- 
However, the structure of the rotation matrix (|47j) implies that 



2k 2l + sin 2 /30 lk O u ~ sin 2 (35 lk , I, k e 2, 3. (48) 

On the other hand, the factor T>ki appearing in the expression for [C3[12]] is nonzero only for 
k ^ I. These two requirements are incompatible, i.e. the total rate asymmetry is strongly 
suppressed.|| 

We thus have to make use of asymmetries that are sensitive to the beam polarization. To 
this end we introduce rate and polarization asymmetries with respect to the unpolarized vector 
boson exchange and Higgs boson exchange parts, {Ci[ij]} and {Csfr/']}, respectively, as follows: 



A P (C 2 [ij]) = Sh 
A R (C 2 [12]) = 2 

a p (c A [ij]) = n 4 



j d cose {c 2 [ij}} 
j dcose{Ci[ijj} ' 

J rfcose[C 2 [12]] 
J rfcos6{Ci[12]} ' 

J d cos Q{C,[ij}} 

/dcose{c 3 [ij]}' 



a (c nan - n LA^1^3 im 

A R (Q112J) - n 4 Jdcos(d{C3m ■ (49) 

The remaining 12 x 4 correlations (n = 5 — 16) can only be measured if at least one of the 
initial beams is transversely polarized. In this case the proper normalization involves the sum 
of the unpolarized Higgs and gauge boson exchange contributions: 



" j dcose ({am + {am) 



^This also implies that in the MSSM the lightest Higgs boson mass must be smaller than the mass of the 
lightest tau slepton pair, so that its contribution to any f pair production process is negligible. 
" This corrects an erroneous numerical result in Ref . |17| . 
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A (C\U])-n I dcosQ[C n [12]] 

A R (C n [12)) - n n j dcosQ{{Cim + {C3m) ■ (50) 



The Q n in Eqs. ( fSf ) and ( |5U| ) are numerical factors that originate from projecting out the corre- 
sponding observable, by adjusting the (anti-)muon polarization and integrating the distribution 
£jj over the azimuthal angle $ with an appropriate weight function, as described below. 

These factors fl n (n = 2, 4, 5 — 16) depend sensitively on the degrees of longitudinal and 
transverse polarizations of the muon and anti-muon beams achievable at muon colliders. For the 
sake of discussion, let us assume that both the muon and anti-muon beams are polarized with 
perfect degrees of longitudinal or transverse polarization. From Eq. (^) we can then deduce the 
algorithms for extracting our asymmetries, which determine the Q n . In the following discussion 
we always assume that all cross sections are appropriately (anti-) symmetrized in the f indices. 

• The asymmetries oc C^fij] can be measured by dividing the difference of cross sections 
for Pl = —Pl = 1 and Pi = —Pl = —1 by the sum of these cross sections; this gives 
Q2 = 1.0 Note that a CP-transformation leaves the polarization states used in this 
prescription invariant. The asymmetries oc Cjij] can be obtained in the same way from 
the cross sections for P L = P L = 1 and Pl — Pl — — 1, so that f2 4 = 1 as well. In this 
case a CP transformation maps the two polarization states into each other; this difference 
of cross sections is therefore a CP-odd quantity, in accord with the results of Table 1. 

• The asymmetries involving C15 and Ciq can be obtained by fixing Pt = Pt = 1 and 
taking the ratio of the difference and the sum of the distributions for cos(a — a) = ±1 
(sin(a — a) = ±1), respectively. (Recall that Pt = Pt = 1 implies Pl = Pl = 0, and that 
a — at is independent of $.) This gives 15 = Qiq = 1. Recall that a CP transformation 
flips the sign of a — a. Since the cosine is invariant under the change of sign of its 
argument while the sine changes its sign, the prescription for extracting C15 is CP-even 
while that for C\q is CP-odd, again in accordance with the CP transformation properties 
of the corresponding terms in the squared matrix element; see Table 1. 

• The asymmetries involving C13 and C14 (which vanishes in the MSSM in the limit V z = 0) 
can also be obtained by chosing Pt — Pt — 1- However, the angular combinations 
cos(a + a) (sin(a + a)) depend on the azimuthal angle see Eq. (p6|). Therefore we 
need to apply a (CP-even) normalized weight function v^2cos(a + a) (v^2sin(a + a)) in 
order to extract these asymmetries, where a + a = rj — 2$. The normalization factors 
in these weight functions /($) are determined by the requirement that / acts as 
a projector, which implies f w d&f 2 (&) = 1. Note that this weight factor should not be 
applied in the denominator of Eq. ( [50]) . Finally, one has to take the sum of these weighted 
integrated cross sections for rj = and 77 = n, in order to remove terms oc C15 or C\q from 
the denominator^ note that these two configurations of transverse polarization vectors 
are self-conjugate under a CP transformation, since r\ = —n is identical to rj = +n. We 



"Recall that [C 2 [12]] = in the MSSM if the Z width is neglected. 

t^The numerator is independent of 77 after the weighted integration over <E>. 
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thus find O13 = fii4 = 1/ v2. Since a + a is also invariant under CP-transformations this 
extraction prescription is CP-even, in accord with the corresponding entries in Table 1. 

• The asymmetries oc C n (n = 5 — 8) are obtained from the cross sections for Pt = 
1, Pt = Pl = and Pt = 1, Pt = Pl = 0, taking r\ = 0. C5 and C$ can be obtained 
from the (CP-even) sum of the corresponding cross sections, while CV and Cs can be 
determined from the (CP-odd) difference of cross sections. In addition one needs the 
weight function \/2cosa; for C 5 and CV, and a/2 sin a for C 6 and Cs] for 77 = these 
weight functions are all CP-even, so that the extraction procedure again has the same 
CP property as the corresponding coefficient of the squared matrix element. This yields 

n 5 = n 6 = n 7 = n 8 = 1/V2. 

• The algorithm for extracting the asymmetries oc C n (n = 9 — 12) is complicated. One 
possible procedure is to measure the distributions for the polarization combinations {Pl = 
±1,Pt = 1} and {Pl = ±l,Py = 1} for rj = 0. Denote these four distributions by E± 
and E± . A CP-transformation then sends E± to E^ and vice versa. The asymmetries 
oc Cg and Cn can then be extracted from E+ — E_ ± (E^ — S^ 2 ^, using the weight 

function \/2cosa; note that the plus (minus) sign in front of the round parentheses gives 
a CP-odd (CP-even) combination of cross sections. The asymmetries oc C\o and oc C\2 
can be extracted in the same way, if the weight function a/2 sin a is used. This gives 
fig = fiio = fin = fil2 = l/v2- 

In reality the degrees of longitudinal and transverse polarization are not perfect, in which 
case the polarization factors fi n should be multiplied by the relevant (products of) degrees of 
longitudinal and transverse polarizations, \ L aric ^ X T (x L an d X T )i °f ^ ne incident muon (anti- 
union) beam. As already stated, we will simply take \ L = X T = X L = X T = ^- m ^ ne following 
numerical analysis. 



5 Numerical results 

In this section, we present some numerical analyses of CP violation both in the tau-slepton and 
neutral Higgs boson systems based on a specific scenario for the relevant SUSY parameters. 
The loop-induced CP violation in the Higgs sector can only be large if both |//| and \A t \ 



(or \Ab\, if tan/? ^> 1) are sizable [|T^, [13], [14|, [15]. For definiteness we will present results for 
tan (3 = 10 and the cm. energy y/s near the two heavy Higgs boson resonances. Since even for 
tan/3 = 10 the contributions from the bottom (s)quark sector are still quite small, our results 
are not sensitive to and A^, we therefore fix Aj, = A t and = = ttLq, although 
different values for the SU(2) doublet and singlet soft breaking squark masses, ttIq ^ m^y, 
are allowed, and also take equal phases for A t and Aj,. Since we are basically interested in 
distinguishing the CP non-invariant Higgs sector from the CP invariant one, we take for the 
re-phasing-invariant phase of A t> &/^e^ two values; (CP invariant case), and 7r/2 (almost 
maximally CP violating case) while the phase of the gluino mass $^ is set to be zero. We chose 
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the following "standard set" of real mass parameters and couplings of the Higgs and squark 
sectors: 

m A = 0.5 TeV, \A t , b \ = 1.0 TeV, = 2.0 TeV, 

m- g = 0.5 TeV, <5> g = 0, nig = m = m b = 1.0 TeV . (51) 

In our numerical analysis, the top and bottom quark running masses m t (m t ) = 165 GeV and 
tfib(mb) = 4.2 GeV will be taken. For rriA much larger than mz as in Eq. (pT|), two neutral 
Higgs bosons have almost degenerate masses: 

= : m H2 = 499.6 GeV, m Hz = 500.0 GeV, 

= - : m H2 = 499.5 GeV, m Ha = 500.4 GeV. (52) 

In the CP invariant case ($A/x = 0) we have cos an = 0, i.e. the state if 2 is CP even and the state 
H 3 is CP odd. On the other hand, our parameter set with &Afi = tt/2 gives cosa# = 0.976, i.e. 
now the state H2 is mostly (but not entirely) CP-odd, and H% is mostly CP-even. Because of 
the near-degeneracy of the two heavy Higgs bosons a significant overlap effect, i.e. interference 
between these two Higgs boson resonances is expected in the presence of CP violation. Note 
that CP violation in the Higgs sector tends to increase the mass splitting between the two 
heavy states; however, in our case it still does not exceed the widths of the two Higgs bosons, 
which amount to about 1.2 GeV. 

We work in the general MSSM, i.e. we do not assume unification of scalar soft breaking 
parameters (masses and A— terms) .[] The real SUSY parameters associated with the tau slepton 
system are thus independent of the parameters in the Higgs and squark sectors. We take as 
standard inputs: 

m fL = 0.23 TeV , m fR = 0.18 TeV , \A T \ = 0.5 TeV . (53) 

The phase <&a t of the rephasing invariant quantity A T [ie 1 ^ is varied. For reference, we list the 
tau slepton masses for $ j4r = 0, n/2: 

$ At = : m fl = 82.7 GeV, m f2 = 287 A GeV, 

$ At = I : m fl = 87.6 GeV, m f2 = 286.0 GeV. (54) 

Since the cm. energy is assumed to be around the resonances of the heavy Higgs bosons, 
y/s ~ 500 GeV, the tau slepton pairs r^ff and f x can be produced in collisions, but 

heavy (f^^ - ) pairs are not accessible, as can be seen from Eq. (^2|). 

In the previous section, we have listed 64 polarization observables which can be constructed 
if the production plane can be reconstructed (at least statistically), i.e. if both the polar angle G 
and the azimuthal angle $ can be measured. We have found that there are 10 CP-violating rate 

^This implies that we cannot make statements about first or second generation sfermion masses. Moreover, 
we do not assume that the electroweak gauge symmetry is broken radiatively. 
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asymmetries and 18 CP-violating polarization asymmetries. For some polarization observables, 
and for all rate asymmetries, it is important to experimentally identify the electric charges of 
the produced sfermions from their decay products. These requirements are very difficult to 
satisfy in the production of third generation squarks. Moreover, in most SUSY models squarks 
are significantly heavier than sleptons, so that squark pair production may not be possible 
for y/s ~ m,A- This is why we focus on the production of tau slepton pairs. Even here the 
requirement that the heavy Higgs boson masses are larger than the sum of the tau slepton 
masses imposes a nontrivial constraint on the accessible parameter space. Moreover, in the 
absence of a detailed Monte Carlo study it is not clear how well the production plane can be 
determined experimentally. In the following we therefore consider the cases without and with 
direct reconstruction of the production plane separately. 

Before turning to these asymmetries we briefly discuss the behavior of the total f+ff" and 
fi?2 cross sections. These total cross sections are not only necessary to estimate the number 
of available events for a given integrated luminosity. Since (apart from overall prefactors) 
they also appear in the denominators of Eqs. ([49|) and (|50|), they also significantly affect the 
•^/i— dependence of the asymmetries. These total cross sections are shown in Fig. || We see 
that the f\ pair production cross section is sizable (~ 100 fb) even away from the Higgs poles. 
Since here the dominant contributions come from gauge boson exchange, there is almost no 
sensitivity to the CP-violating phase $>a t -> as discussed earlier. On the other hand, the Higgs 
exchange contribution to the total cross section does depend significantly on this phase. Not 
only the maximal value of this cross section, but also the energy where this maximum is reached 
depends on the CP-violating phases. This follows from the observation that in the absence of 
CP-violation a CP-odd Higgs boson cannot couple to an identical f pair; however, once CP is 
violated, both heavy Higgs bosons contribute. In fact, for &a t — tt and y/s ~ m^, Higgs boson 
exchange clearly dominates the total f\ pair production cross section. 

The dominance of the Higgs boson exchange contributions near the poles is even more 
pronounced in the f\fi channel.0 Our choice (|53"D implies |sin2# T | ~ 1, which maximizes 
the Zf\fi coupling fll8|) . Nevertheless the gauge boson exchange contributions are suppressed 
compared to the case of f\ pair production by the absence of photon exchange contributions 
and by the smaller available phase space (/3f 2 — 0.24, as compared to (5\ x ~ 0.83). Note that 
the Higgs boson exchange contribution only scales like f3, and is hence much less phase space 
suppressed. Moreover, even in the absence of CP-violation both Higgs bosons can contribute 
to f\fi production. Nevertheless a strong dependence on the value of can be seen also in 
this case. Note that the fif 2 production cross section typically changes by more than a factor 
of 10 over the shown range of energies, as compared to a variation by a factor 2 to 5 in case 
of fi pair production. The energy dependence of the total cross section is therefore especially 
important for the understanding of the a/s dependence of asymmetries in the f\fi channel, 
which includes all rate asymmetries. Note also that the peak cross section, which falls in the 
range of ~ 50 to ~ 300 fb, is sizable also in this channel, when compared to the anticipated 
luminosity of ~ 10 fb -1 per year [22|. Finally, we note that in the absence of CP-violation 



^The difference between the totai cross sections for f^f 2 and f^f x production is negligible, as shown in 
Sec. 4.2. 
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in the Higgs sector the total cross sections depend are even functions of the phase $a t - This 
explains why the curves for <&a t = <fi coincide with those for = 27r — in Figs. 0(a) and 
(c). However, this degeneracy is lifted if 7^ 0. 
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Figure 3: The yfs dependence of the total fx pair cross section (a,b) and the mixed t+tvT cross 
section (c,d) near the heavy Higgs boson resonances, with the SUSY parameter set of Eqs. (pjl), 
( [5% ) and (\5$). Frames (a) and (c) are for = 0, while frames (b) and (d) are for $ a^ = 7r/2, 
The six lines in each frame are for $a t = (solid), n/3 (long dashed), 2n/3 (dot-dashed), n 
(dotted), Att/3 (short dashed) andhn/3 (heavily dotted), respectively. 



5.1 Without direct reconstruction of the production plane 

If the production plane is not reconstructed, all the SV correlations are averaged away when 
integrating over the azimuthal angle $, but the polarization and rate asymmetries corresponding 
to the coefficients {Cx, C2; C3, C4, Cxs, Cxe} survive. As emphasized earlier, the VV correlations 
are, however, determined by the SUSY parameters associated with the tau slepton system 
without any CP-preserving phases, if the Z boson width is neglected. In this case the CPT- 
odd VV rate asymmetries ^t^(Ci j 2[12]) vanish so that there exist no non-trivial CP-violating 
VV asymmetries^ for probing the CP phases in the tau slepton mass matrix. This immediately 

*From now on, let us denote the asymmetries with SS, VV and SV correlations in the numerators by SS, 
VV and SV asymmetries, respectively. 



25 



implies that f pair production in e + e~ collisions is of no use for directly measuring the CP 
phase associated with the tau slepton system.^ 

On the other hand, the SS polarization and rate asymmetries are nonzero. The 6 CP- 
violating SS polarization asymmetries Ap(Ci[ij]) and Ap(Ci 6 [ij}) depend on CP violation in 
the f sector only through the CP-even combination of couplings ^(V^V^) , see Eq. (f44|); 
these asymmetries are CP-odd, since they are also proportional to the CP-odd combinations 
Oik02iT02kOu of Higgs mixing angles. The product of this combination of Higgs— ff couplings 
and Higgs mixing angles vanishes if CP is conserved, but in general one can expect a nonzero 
result in the presence of CP violation in either the Higgs or f sector. If there is no CP- 
violation in the Higgs sector = 0), these (and all other CP-violating) asymmetries are 
odd functions of sin$ J 4 T , i.e. the transformation —* 27r — leads to a change of sign 
of these asymmetries, leaving their absolute values unchanged. Since these are polarization 
asymmetries, they should be nonzero even for the production of two equal f states, in particular 
for fi pair production. 

These expectations are borne out by the numerical results shown in Fig. Here the upper 
two frames show Ap(C4[ll]) and the lower two frames show t4p(Ci6[11]), while the left (right) 
frames are for a CP-conserving Higgs sector ($a^ = tt/2). Since Ap(C4[ll]) is CPT-odd, 
it is proportional to the imaginary part of the product of relevant Higgs boson propagators, 
while A p (Ciq), being CPT-even, is proportional to the real part of the product of propagators. 
Eq. ( f£2"D then explains why y4p(C4[ll]) decreases much faster as one moves away from the 
Higgs poles than v4p(Ci6[ll]) does. The results of Fig. [| show that both of these polarization 
asymmetries are quite sensitive to A T . Also, both asymmetries can be large. Finally, we note 
the curious fact that in the presence of sizable CP-violation in both the Higgs and f sectors 
some CP-odd asymmetries can be "accidentally" suppressed; see the long dashed curve (for 
$ AT =vr/3)inFig.|(b). 

The possibly large size of Ap(C4[ll]) is very encouraging, since this asymmetry can be mea- 
sured with only longitudinally polarized beams; this is important, since it might be technically 
quite difficult to produce transversely polarized muon beams. Moreover, this asymmetry is lin- 
ear in polarization; this means that the precision with which it can be measured also depends 
"only" linearly on the achievable degree of polarization. Moreover, in the presence of large 
CP-violating phases this asymmetry remains sizable over a large region of parameter space. 
This is shown in Fig. |5|, which shows the behaviour of the maximal value of Ap(C4[H]) when 
one input parameter is moved away from our standard input set described at the beginning of 
this Section. The maximum here refers to the dependence on i.e. we plot the asymmetry at 
the optimal center-of-mass energy, which depends (slightly) on the input parameters. This is 
legitimate, since i/i can be freely chosen in an experiment. In order to show the dependence on 
different parameters within a single frame, we chose as x— axis the ratio of a given parameter to 
its "standard" value. The range of variation of these parameters is restricted by experimental 
constraints, the most important ones being the lower bounds rriH 1 > HI GeV (note that Hi 
behaves essentially like the Standard Model Higgs boson, since m 2 A 3> Mf ) and m ?1 > 80 GeV. 

t An exception to this rule can occur if the two f mass eigenstates happen to be nearly degenerate, with mass 
splitting comparable to their decay widths [E9| . 
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Figure 4: TTie -^i dependence of the CP-odd SS polarization asymmetry Ap(C4[ll}) for (a) 
^Afi = and (b) §Afi = tt/2, and the CP-odd SS polarization asymmetry ^4p(Ci6[ll]) for (c) 
&Afi = and (d) = 7r/2. Parameters and notation are as in Fig. 3. 



We see that the maximal asymmetry depends only rather weakly on tan/? (dashed curve). 
However, the total cross section near the Higgs resonance falls oc tan 4 (3 for tan 2 (3 m t /mb, due 
to the decrease of the couplings of the heavy Higgs bosons to f sleptons and the simultaneous 
increase of the Higgs decay widths into ti quarks. The dependence on \A T \ is also very weak, as 
long as \A T \ ^200 GeV. However, sending \A T \ — > removes CP-violation from the f sector and 
also greatly reduces the couplings of Higgs bosons to f pairs, see Eq. (|20|) ; hence all CP-odd 
asymmetries become very small in this limit. The dependence on \fi\ (dot-dashed curve) is 
somewhat stronger than that on \A T \, however, we again find a maximal asymmetry £ 0.5 
unless \n\ is reduced by more than a factor of 3 from its standard input value of 2 TeV. Finally, 
the maximal asymmetry remains finite even as \A t \ = \Ab\ — >• 0, as shown by the dotted curve. 
In this limit there is no CP-violation in the Higgs sector, but we already saw in Fig. (| that 
CP-violation in the f sector by itself is sufficient to produce a large value of this asymmetry. 
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Figure 5: The maximal value of \Ap{Ci[ll})\ is shown for = tt/2 and §a t = 47r/3 7 where 
"maximal" means that the optimal value of y^s has been taken. The other input parameters 
have been fixed to their "standard" values, except that one of these parameters has been varied 
along each curve. The x—axis is the ratio of this parameter to its "standard" value: for the 
solid curve, R = \A T \/(0.5 TeV); for the dashed curve, R = tan/3/10; for the dot-dashed curve, 
R = |/i|/(2TeV); and for the dotted curve, R = \A t \/(l TeV). The upper bounds on tan/3 and 
come from the requirement m fl > 80 GeV, while the lower bound on tan/3 comes from the 
bound m Hl > 111 GeV. 

Unfortunately studies of f\ pair production alone do not suffice to completely determine 
the parameters of the f sector. Even under the most favorable circumstances such studies can 
only determine the values of four real parameters: the mass , the mixing angle 6 T (these 
two parameters can already be determined from VV correlations), and the couplings of the 
two heavy Higgs bosons to fx pairs. On the other hand, the f mass matrix depends on 6 real 
parameters (m fL ,mf R , \A T \, |//|,tan/3 and $a t )- Moreover, if the mixing between CP-even and 
CP-odd Higgs states is small, C 4 [ll] and Ci 6 [ll] essentially only depend on the product of the 
couplings of the two heavy Higgs bosons to fi pairs, since in such a situation the product of 
mixing angles 0\kOi\ will be sizable only for k ^ I. In our numerical example this is trivially 
true for = 0, but is also approximately correct for $^ = 7r/2, where H 2 is predominantly 
CP-odd while H 3 is mostly CP-even. Fortunately the total cross section for fi pair production 
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for y/s ~ vriA should depend on a different combination of Higgs couplings, so that it might 
still be possible to determine the couplings of both heavy Higgs states to f\ pairs. 

If mixed f^f^ pairs are accessible, many more quantities become measurable. We saw above 
that studies of f\ pair production should allow to determine 4 real parameters in the f sector. 
In principle the total f^f^ cross section, measured in the vicinity of the heavy Higgs poles, 
should then be sufficient to completely fix the parameters of the f sector. However, this is only 
true if we assume that the masses and mixing angles of the neutral Higgs bosons are known. 
In addition, we have to assume that f pair production is indeed described by the MSSM, and 
that there are no other diagrams contributing; for example, in the presence of slepton flavor 
mixing, t— channel diagrams might contribute [p4]] . Analyses of the new asymmetries that 
become available if f 1 f2 pairs can be produced can then be used to check the consistency of 
the framework followed in this paper. Moreover, (some of) these new asymmetries should also 
have quite different systematic uncertainties than the quantities discussed so far. 

In particular, we now have the opportunity to study rate asymmetries in addition to polariza- 
tion asymmetries. All the SS rate asymmetries, i.e. the two CP-even asymmetries ^4^(64 [12]) 
and v4r(Ci 6 [12]), and the two CP-odd asymmetries ^t^(C 3 [12]) and -4r(Ci 5 [12]), are in prin- 
ciple sensitive to CP violation in the tau slepton system even without reconstruction of the 
production plane. However, we already saw at the beginning of Sec. 4.2 that [C*3[12]] is always 
very small in the MSSM with m\ 3> m 2 z . Moreover, we find numerically that v4r(Ci6[12]) is 
relatively insensitive to <&a t , m particular if there is no CP-violation in the Higgs sector. In 
Fig. I we therefore show A R {Ci [12]) (top row) and t4#(Ci 5 [12]) (bottom row). 

The first of these asymmetries is even under both CP and CPT. It is therefore nonzero 
even if = = 0, and also remains finite far away from the heavy Higgs poles; recall that 
this asymmetry is normalized to the (unpolarized) squared Higgs boson exchange contribution 
only, see Eq. ([49]). However, this is only true for completely (transversely) polarized beams; 
and even in this case the number of events will become very small as one moves away from 
the Higgs poles. Nevertheless measuring this asymmetry in the vicinity of the pole should 
yield useful information. Note that in the absence of CP-violation in the Higgs sector, the 
asymmetry depends only on cos<&a t , i.e. it becomes invariant under — ► 2n — &a t ', this is 
true for all other CP-conserving asymmetries as well. However, this degeneracy is broken if 
7^ 0. Generally we find that CP-violation in the Higgs sector increases the sensitivity to 
the rephasing invariant phase in the f sector. 

This is true also for y4^(Ci5[12]). In fact, this asymmetry is very small if there is no CP- 
violation in the Higgs sector. The reason is that the products of couplings C^fcC^ and OikOu 
appearing in this asymmetry, see Eq. ([4"4|), are almost proportional to Ski if the heavy CP-even 
and CP-odd Higgs states do not mix, while the combination TD^i of Higgs boson propagators 
vanishes for k = I. The interference between the heavy and light CP-even Higgs states then 
at best gives an asymmetry of order 10~ 3 , which is not measurable with realistic luminosities. 
On the other hand, once CP is violated in the Higgs sector, this asymmetry can become large. 
However, since it is CPT— odd, it decreases quickly as we move away from the Higgs poles. 
Moreover, the results of Fig. |6] indicate that this asymmetry is not very sensitive to even 
if'1% / 0. 
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Figure 6: The y/s dependence of the CP-even SS rate asymmetry ^4^(C 4 [12]) for (a) $^ = 
and (b) $^ = ir/2, and the CP-odd SS rate asymmetry Ar(Ci 5 [12}) for (c) $a = and (d) 
&a — 7r/2. Parameters and notations are as in Fig. 3. 



5.2 With direct reconstruction of the production plane 

If the production plane of the tau slepton pair can be reconstructed efficiently, the SS asym- 
metries related to C13 and C14 become accessible. More importantly, all SV correlations can 
now be studied. These asymmetries are linear in Higgs propagators, which implies that they 
are less sensitive to the strength of the Yukawa coupling of the muon than the SS correlations 
are. Moreover, in principle asymmetries oc Kef/^J allow one to "switch off" the contribution 
from Higgs mass eigenstate Hk simply by setting ^/i = m# fc , since the real part of the Higgs 
propagator vanishes there. In this manner one might be able to cleanly isolate the contribution 
from the second heavy Higgs boson. In this section, based on the parameter sets (|5T|) and 
(|5BJ), we investigate quantitatively whether the SV asymmetries can provide us with additional 
useful information on the f sector. 

In principle SV correlations can already be studied in polarization asymmetries in the rj^rf 
channel. Unfortunately we find numerically that these asymmetries are always rather small, 
with absolute value ^ 0.2. This is true also for polarization asymmetries in the mixed f^f^ 
channel. We are thus left with the 8 rate asymmetries ^4^(C n [12]) (n = 5 - 12). Eq. ( ^3|) shows 
that all these asymmetries depend on the parameters of the f sector through the combination 



30 



of couplings 



h = 3m (V k . 12 Qf*) , (55) 

as well as through the masses of the f sleptons. A closer look at the explicit expressions for 
these rate asymmetries shows that the four asymmetries Af?(Cn[12]) (n = 6, 7, 10, 11) are 
proportional to the small numerical factor 4s^/ — 1 ~ —0.07 for ~ 0.233; this comes from 
the vector coupling of the Z boson to muons. We will therefore not consider these asymmetries, 
and instead concentrate on the remaining four asymmetries, which stem from the axial vector 
coupling of the Z boson to muons. They satisfy the proportionality relations: 



A R {C 5 [12]) <x-0 2k $m[D Hk ]I k , 
A R (C 8 [12]) <x+sp0 lk Sm[D Hk )I k 
A R (C g [12}) cx +sp lk fte[D Hk ] h, 
A R (C 12 [12])<x-0 2k $te[D Hk \I k . 



(56) 
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Figure 7: The \fs dependence of the CP-violating SV rate asymmetry .4r(C 5 [12]) for (a) 
^a^i = and (b) = n/2, and the CP-conserving SV rate asymmetry Ar(C 8 [12]) for (c) 
^a^i = and (d) $^ = tc/2, respectively. Notations and parameters are as in Fig. [|. 



These relations show that the CP-conserving asymmetries Ar(Cs,9) involve the CP-odd com- 
ponents (Oik) of the Higgs bosons, while the CP-violating asymmetries Ar(C^i 2 ) involve the 
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CP-even components (02k) of the Higgs bosons. This means that in the absence of CP violation 
in the Higgs sector, only one of the two heavy Higgs bosons contributes to a given asymmetry. 
In our case this remains approximately true even if 7^ 0, since, as noted earlier, the quan- 
tity sin(2ct#), which measures the strength of CP-violation in the Higgs sector, is quite small. 
In other words, in our examples one does not even have to tune y/s ~ m H k m order to isolate 
the contribution of a specific Higgs boson. 

Note also that the two CPT-odd asymmetries Ar(C5,8[12]) have their peaks at the heavy 
Higgs boson poles, but they are suppressed far away from those poles due to 3m[D^J. However, 
Fig. [7| shows that these asymmetries only decrease relatively slowly as one moves away from the 
poles. This can be understood from the strong energy dependence of the total f\Ti production 
cross section shown in Fig. |3|: the denominator of ^4r(C 5i8 [12]) also decreases quickly as one 
moves away from the poles. 




Figure 8: The y/s dependence of the CP-conserving SV rate asymmetry Ar{C$[\2}) for (a) 
^Afi = and (b) $,4^ = 7r/2, and the CP-violating SS rate asymmetry Ar(Ci2[12]) for (c) 
^Api = and (d) = tt/2, respectively. Parameters and notation are as in Fig. [|. 



On the contrary, Fig. H shows that the other two asymmetries •Ai?(Cg i i2[12]), which are 
CPT-even, reach their extrema several GeV away from the Higgs poles. The numerators of 
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these asymmetries actually reach their extrema at ^/s = niu ± r#/2, where m# and Yh are 
a typical heavy Higgs mass and decay width, respectively. However, the rapid decrease of 
the denominator pushes the extrema further away from the Higgs poles. Note also that these 
asymmetries can attain larger values than the CPT-odd asymmetries shown in Fig. |7|. How- 
ever, while .Ar(C5 ) 8[12]) can be measured if only one initial beam is (transversely) polarized, 
Ar(C9,i2 [12]) can only be determined if both muon beams are polarized, one having a trans- 
verse and the other a longitudinal polarization. Together with the requirement that the event 
plane has to be reconstructed this means that measuring the asymmetries shown in Fig. |] will 
probably pose the biggest challenge to both collider and detector. 



6 Summary and Conclusions 

We have performed a detailed, systematic investigation of the signals for CP violation in the 
neutral Higgs boson and tau-slepton systems through the production of tau slepton pairs in 
polarized collisions, — > f^fj~ with the labels i,j = 1,2 for the two r slepton mass 

eigenstates. We worked in the framework of the MSSM with exact R-parity and negligible 
flavor mixing. The relevant sources of CP violation can then be found in soft breaking terms 
associated with third generation sfermions, as well as the \l— parameter. CP violation in the 
f sector contributes to CP-odd asymmetries at the tree-level. CP violation in the t and/or b 
sector leads to mixing between CP-even and CP-odd Higgs current eigenstates. Even though 
this CP-violating Higgs mixing only proceeds through loop diagrams, it can give rise to 0(1) 
CP-violating asymmetries even in the absence of other CP phases. 

The expression for the cross section for f pair production from a yU + yU~ initial state with 
arbitrary (possibly transverse) polarization contains 16 terms. We have classified the behavior 
of these terms under CP and CPT transformations. Terms with CP-even polarization factors 
can contribute to CP-violating rate asymmetries, i.e. differences between the cross sections 
for fi~T2 and fff^ production. On the other hand, terms with CP-odd polarization factors 
contribute to CP-violating polarization asymmetries that can be probed already in fi pair 
production. In some cases it can also be interesting to construct asymmetries from these terms 
which are even under a CP transformation. These asymmetries are nonzero even in the absence 
of CP violation, but they can nevertheless help to determine the to date unknown parameters 
in the problem, including phases. 

It is reasonable to assume that all the properties of the neutral Higgs bosons are determined 
beforehand, e.g. through the study of n + fi~ — > ff production near the Higgs poles, where 
/ = r,b or t. The remaining task is thus the determination of the parameters in the f sector. 
As a first step one might want to determine the masses of the accessible f states, as well as the 
?l — tr mixing angle 9 T through the study of f pair production away from the Higgs poles. Here 
the cross section will be dominated by gauge interactions, and is thus completely determined by 
these parameters. The mixing angle 9 T also affects the polarization of the r leptons produced 
in f decays J3"5|. Of course, these measurements can also be performed at an e + e~ collider. 
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Figure 9: The yfs dependence of the integrated luminosity required to detect the asymmetries 
Ap(C4[H]) ( a ,b) and Ar(C 5 [12]) (c,d) with statistical significance of one standard deviation, 
for ideal beam polarization. The luminosity increases quadratically with the required number of 
standard deviations, and with the inverse of the degree of beam polarization. Parameters and 
notation are as in Fig. [|. 



f pair production through gauge interactions can thus at best give three relations among 
the six real parameters appearing in the f mass matrix. The remaining parameters can be 
determined if there are sizable Higgs exchange contributions. In particular, the longitudinal 
polarization asymmetry ^4p(C 4 [ll]) is usually ideally suited for directly probing CP-violation 
in the f sector. One does not need to reconstruct the event plane to measure this asymmetry, 
nor are transversely polarized beams required. Moreover, this asymmetry is linear in the 
achievable polarization; this is of some importance, since the actual beam polarization might 
be significantly less than the ideal value of 100% assumed in the numerical results of Sec. 5. 
Of course, we have to assume that f pair production is kinematically allowed in the vicinity of 
the Higgs poles. Moreover, the ratio of vevs tan (5 must not be too small, since the couplings of 
both muons and tau (s)leptons to Higgs bosons scale essentially oc tan/3. However, scenarios 
with very small tan (3 are essentially excluded by LEP Higgs searches anyway. The upper two 
frames in Fig. |9] show that, if these conditions are met and the CP-violating phases are not 
very small, this asymmetry should indeed be measurable with the foreseen luminosity of about 
10 fb _1 per year. Here we show the integrated luminosity required to see a lex deviation of this 
asymmetry from zero, for ideal beam polarization. The required luminosity scales quadratically 
with the required number of standard deviations as well as with the available polarization. For 
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example, a value of 0.01 fb _1 means that with a data set of 10 fb _1 and 35% beam polarization 
one can see an Ilex effect, i.e. this asymmetry would be measurable with relative precision of 
about 9%. We saw in Fig. [5] that this conclusion should remain true over a wide region of 
parameter space. 

Note, however, that this asymmetry might be "accidentally" suppressed, even though there 
are large CP-violating phases in both the f and t sectors. It is therefore important to try 
and measure several independent CP-violating asymmetries. Moreover, analyses of f\ pair 
production in the MSSM can at best determine the couplings of f\ pairs to the two heavy 
neutral Higgs bosons. LEP search limits on together with the bound mu 1 ~ 130 GeV 
imply that the exchange of the lightest neutral Higgs boson H\ contributes negligibly to f 
pair production. Simple parameter counting then implies that, in the absence of any prior 
knowledge of the parameters of the mass matrix, one will have to study T\fi production in 
order to completely determine the parameters of the f mass matrix. Note that some of these 
parameters (tan/3 as well as /i) also appear in other terms in the MSSM Lagrangian. It is 
thus possible that their values will be known (with some error) beforehand. Even in that case 
it would be important to determine them independently from f pair production, in order to 
confirm that the MSSM with the given assumptions can indeed describe this process. Of course, 
this will only be possible if rriA > rrif 1 + . 

Unfortunately we found that CP-violating asymmetries that can be measured in the mixed 
f\f% channel are smaller than *Ap(C4[ll]), at least in the vicinity of the Higgs poles where 
the event rate is sizable. This is demonstrated by the lower two frames of Fig. [|, which 
shows the luminosity required to see a la deviation in the rate asymmetry .Ar(C5[12]). In 
most cases the minimal required luminosity is at least an order of magnitude larger than for 
ApCCJH])- We saw in Sec. 5 that CPT— even asymmetries can be 0(1) even in the mixed f\fi 
channel. However, these asymmetries are proportional to the real parts of heavy Higgs boson 
propagators, which vanish near the poles. As a result, the minimal luminosity required for 
measuring these asymmetries (not shown) is similar to that required to measure v 4k(C < 5[12]), 
but would have to be taken at ^fs ~ m#±r#, where m# and Y n are the typical mass and total 
decay width of the heavy neutral Higgs bosons. It thus seems unlikely that studies of mixed 
f\f% production can contribute significantly to the determination of the CP-violating phase in 
the f mass matrix. Fortunately the determination of all parameters of this mass matrix can 
be completed by measuring a single quantity in the f\fi channel, for example the total cross 
section near the Higgs poles. We saw in Sec. 5 that this cross section is indeed often sizable, and 
should thus be readily measurable. A possible problem here might be the distinction between 
f\ pair and fif 2 production. This should be fairly easy in the scenario we studied in Sec. 5, 
due to the large f 2 — f x mass splitting, but could be more problematic if this mass splitting is 
small. 

We thus see that in the chosen framework transverse beam polarization is not necessary 
to determine all the free parameters of the f system. However, in order to test the model 
one has to over-constrain it, i.e. there should be more measurements than parameters. Also, 
generalizations of the model, e.g. allowing R-parity violation and/or slepton flavor mixing, are 
conceivable. In that case the squared matrix element can still be written as in Eq. (|3"3"D, but 
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36| to the expressions 



and hence to the coefficients 



there will be additional contributions 
C n . It is therefore important to measure as many different asymmetries and distributions as 
possible, including those that can only be accessed with transversely polarized beams. Only 
then will it be possible to fully exploit the physics potential of muon colliders to probe the 
details of supersymmetric models through scalar r production. 
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